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Summary 

A turbulence model has been developed for blade- 
element helicopter simulation. This model, called Simu- 
lation of Rotor Blade Element Turbulence (SORBET), 
uses an innovative temporal and geometrical distribution 
algorithm that preserves the statistical characteristics of 
the turbulence spectra over the rotor disc, while providing 
velocity components in real time to each of five blade- 
element stations along each of four blades. 

An initial investigation of SORBET has been performed 
using a piloted, motion-based simulation of the Sikorsky 
UH60A Black Hawk. Although only the vertical compo- 
nent of stochastic turbulence was used in this investiga- 
tion, vertical turbulence components induce vehicle 
responses in all translational and rotational degrees of 
freedom of the helicopter. 

The single-degree-of-freedom configuration of SORBET 
was compared to a conventional full 6-degrees-of- 
freedom baseline configuration, where translational 
velocity inputs are superimposed at the vehicle center of 
gravity, and rotational velocity inputs are created from 
filters that approximate the immersion rate into the 
turbulent field. For high-speed flight the vehicle 
responses were satisfactory for both models. Test pilots 
could not distinguish differences between the baseline 
configuration and SORBET. In low-speed flight the 
baseline configuration received criticism for its high fre- 
quency content, whereas the SORBET model elicited 
favorable pilot opinion. 

For this helicopter, which has fully articulated blades, 
results from SORBET show that vehicle responses to 
turbulent blade-station disturbances are severely attenu- 
ated. This is corroborated by in-flight observation of the 
rotor tip path plane as compared to vehicle responses. 

The Simulation 

A piloted simulation was performed on NASA Ames 
Research Center’s Vertical Motion Simulator (VMS) in 
1994 using the GENHEL UH60A Black Hawk math- 
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ematical model of reference 1 as a platform to investigate 
the Simulation of Rotor Blade Element Turbulence 
(SORBET) model. This was a simulation technology 
experiment designed for the purpose of evaluating the 
influence of stochastic turbulence created at the blade- 
element stations of a rotor system. 

This project was conducted primarily to determine if a 
representative workload could be accommodated by our 
computer system in real time, and to ascertain whether 
the cyclostationary (periodically stationary) effects of 
reference 2 were significant in the motion environment 
of the VMS. Favorable pilot opinion was somewhat sur- 
prising because independent longitudinal and lateral 
inputs were not used during the piloted runs. Furthermore, 
angular components of turbulence were not explicitly 
included in SORBET, but rather occurred as induced 
responses from the superposition of vertical components 
at the blade stations. 

The complete SORBET-BIack Hawk model, including all 
the requisite motion and visual communication software, 
was executed on a single processor (VAX 9000) with a 
cycle time of 12 milliseconds. Cyclic phenomena are 
generally experienced in blade-element rotor models. The 
introduction of turbulence did not amplify these effects. 

For the initial proof-of-concept simulation, only the 
vertical component of turbulence was used. Although 
pilots recognized the absence of the horizontal distur- 
bances under certain flight conditions, their comments 
were generally quite favorable concerning the vehicle 
performance in turbulence. As opposed to the conven- 
tional body-fixed approach to turbulence modeling, 
SORBET has, at least incrementally, improved realism. 

The SORBET Model 

For implementation into the GENHEL UH60A mathe- 
matical model of reference l, an heuristic model of 
turbulence (SORBET) was developed that preserves the 
Gaussian statistical characteristics of turbulence filters 
over the rotor disc. This model avoids the large com- 
putational and storage requirements of recent investiga- 
tions into rotorcraft turbulence models (refs. 2 and 3), and 
imposes a modest computational overhead. SORBET was 
designed to create and propagate all three translational 



components of MIL SPEC turbulence (ref. 4) in real time. 
The turbulence velocity components are superimposed at 
each blade-element location as a function of vehicle 
velocity, with distribution according to the element’s 
instantaneous geometry. 

Recent advances in the computational capacity of the 
simulation facility were significant in the decision to 
develop a model that distributes turbulence over the rotor 
disc. This was anticipated in reference 5: 

Progressive gust penetration by the rotor is an 
important feature which should be included, where 
possible, as it is likely to strongly influence the simu- 
lation results. Consequently, the individual blade 
mode with the gust field model should be used when 
high computation time is permissible. 

The entire stochastic rotor-element turbulence model is 
described here, although for the initial proof-of-concept 
simulation only the vertical component of turbulence was 
used. Examples from this model are presented for the 
helicopter in low-speed flight (10 knots). 

In this study the Dryden spectral form was selected due to 
its computational efficiency, although a discrete model of 
a curve-fit approximation to the von Karman form 1 does 
not offer any particular technical challenge. As will be 
discussed, the rotor-to-body attenuation is so severe that 
the original gust spectra appear to be of secondary con- 
sideration. From reference 4, “When no comparable 
structural analysis is performed or when it is not feasible 
to use the von Karman form, use of the Dryden form will 
be permissible.” The difference between the Dryden and 
von Karman form is basically a small variation in the high 
frequency content (ref. 6). 

SORBET Model Assumptions 

Conventional MIL SPEC angular turbulence filters are a 
substitute for the finite-element distribution of transla- 
tional turbulence in the air mass over the airframe. In the 
current study only the rotor disc is considered as an 
airframe, and the translational turbulence in the air mass 
is distributed over the rotor disc to each blade element. 
This distribution produces rotor moments as well as 
forces. Hence, only the MIL SPEC translational turbu- 
lence filters were used in the SORBET model. 

From Taylor’s hypothesis, the turbulent flow is stationary 
and homogeneous in the longitudinal direction. Turbu- 
lence filter outputs are statistically valid at either 


^he earliest edition of such a model that we have discovered is 
in an unpublished Boeing document (YC-14-FC-1 1 1R) by 
James H. Vincent in December 1973. 


fixed spatial locations or at points translating with the 
vehicle (as is the case in the conventional body-fixed 
formulation). Two points of applicability are selected 
along an onset line that is perpendicular to the vehicle 
flightpath, and tangential to the leading edge of the 
rotor disc. 

The outputs of turbulence filters are applicable at points 
on the onset line at the discrete time points. Interpolation 
perpendicular to the onset line requires a transport delay 
that is proportional to the distance along the flightpath to 
the element location, and inversely proportional to the 
vehicle’s aerodynamic velocity. 

Simulating transport delays requires that the turbulence 
filter outputs be tabled, and updated each computer cycle. 
A practical limit on the table size produces a corre- 
sponding minimum velocity of propagation through the 
rotor disc in simulation. 

The MIL SPEC translational filters contain a pole pro- 
portional to the aerodynamic velocity. This velocity is 
limited to the value determined from the selected 
table size. 

The turbulence is assumed uncorrelated at each side of 
the rotor disc, and homogeneous in the lateral dimension. 
The statistical properties of the turbulence are then 
invariant over the entire rotor disc by using what we call 
Gaussian interpolation between values applicable on each 
side of the rotor disc. 

Geometrical Considerations 

The Black Hawk rotor system has four blades (N = 4). For 
each of the blades the number of segments used in the 
simulation was five ( M = 5). The blade and segment 
indices are defined 

I < n < N (blades) 

1 < m < M (segments) 

Using the equal-annuli algorithm of reference 1, the radii 
to the blade stations may be computed from 



where the rotor radius is R = 26.83 ft, the hinge offset is 
e = 1 .25 ft, and the spar length is e' - 2.25 ft. 

The azimuth angle of the nth blade with respect to the 
vehicle’s aft centerline is given by 
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where Q-21 rad/sec is the nominal rotor RPM. 

In the SORBET model, rather than defining the trans- 
lational disturbance velocities at the vehicle center of 
gravity, they are defined at two onset points, located 
along an onset line that is perpendicular to the horizontal 
mean wind vector v H shown in figure 1 . Using Gaussian 
inputs, three independent turbulence velocity histories are 
computed at each of these two points, in real time. These 
histories are updated each cycle time, and the pertinent 
values from these histories are distributed to the indi- 
vidual segment locations by temporal and geometrical 
algorithms. 



Where u h is the vehicle’s longitudinal aerodynamic 
velocity and v h is the lateral aerodynamic velocity, the 
in-plane aerodynamic velocity may be defined 
(approximately) by 

I 2 , ~2 
V H=i“ b + V b 

Using this velocity and the sideslip angle /J, the perpen- 
dicular onset line is as shown in figure 1 . The aerody- 
namic azimuth angle (wind axis) tor the nth blade is then 

Y n =Vn + P 

where y/ n is the nth blade’s geometrical azimuth angle 
from aft. 

This geometry is shown below to be sufficient to establish 
both the longitudinal and lateral distributions. 


Longitudinal Distribution 

The left and right Dryden filters receive Gaussian inputs, 
so the outputs of these linear filters are also Gaussian. 
Furthermore, a transport delay in these outputs does not 
alter the statistical characteristics. The outputs of the 
filters are stored in tables, which are updated each cycle 
time. The value applicable for a given element is deter- 
mined by selecting a value from the tables. Given that the 
table size is established ( K M = 500) for storing time 
histories of the outputs of the filters, the minimum aero- 
dynamic velocity that may be accommodated is 
determined by distributing the table over the rotor 
diameter. 

2 R 

Vmin = K M A, 

If a smaller minimum velocity is required, a larger table 
size must be selected. For the SORBET simulation cycle 
time A t = 0.012 sec, the minimum velocity for K M = 500 
is 8.944 ft/sec (5.3 knots). The rotor-plane aerodynamic 
velocity used in the turbulence filters is restricted to this 
minimum. 

[ v min V H < v min 

^UV ~ I 

{ V H V H - v min 

Independent of vehicle velocity, exactly K M At = 6 sec of 
data are stored in the tables. Time histories of the outputs 
of six turbulence filters are contained in six separate 
tables, each of length K M . At the minimum velocity the 
individual cells correspond to a length of 0.10732 ft (such 
that 500 of them span the rotor diameter), whereas at a 
velocity of 100 ft/sec (59.2 knots), for example, the cells 
are each 1.2 ft long such that only 45 of them are needed 
to span the rotor diameter. 

In order to select the velocities to be used from the tables, 
for any blade element an integer index is computed to 
determine the element’s perpendicular distance from the 
onset line. Defining 

r m,n = R + r m cos 
the index is given by 

, _ [ r mj j 1 


where the upper brackets denote the integer ceiling 
operation. The tables extend beyond the rotor disc for 
higher vehicle velocities than the minimum, such that 
stationary outputs are available for table size increases 
caused by decreases in velocity. The hypothetical case 
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where the entire rotor disc receives the same table value 
could only occur for velocities in the thousands of knots. 


Lateral Distribution 

Having determined the longitudinal point of application 
(or index k m n ) for a given rotor element, the turbulence 
velocities from the tables are identified as applicable at 
both the left and right sides of the rotor disc. Lateral 
interpolation must be used to determine the element’s 
final turbulence velocities. Although the lateral pro- 
portion (from the left-hand side) to an element may then 
be determined from the geometry of figure 1 

1 . Vn siny »i 
p m.n 2 2R 

this proportion cannot be applied in a linear fashion. 
Using the vertical dimension, for example, although both 
Aw L (k m n ) and Aw R (k m n ) are Gaussian random variables 
with zero mean and RMS value given by C H ., in order to 
preserve the statistical properties over the rotor disc we 
use Gaussian interpolation, so that the combination 
retains the same variance (and all other moments) as the 
original Gaussian variables. 

A variable x is normally distributed with zero mean and 
variance cf 2 when its density function is 

_ r mM s * n 
Pm.n ~ j JR 


P x 1 


0 - p) x 2 


-J/? 2 +(l-p) 2 -Jp 2 +( 1 - p) 2 

P 2 * I +(1 ~ pfx 2 + 2 p(l - p)x { x 2 


p +0 -p) 


p 2 * 2 + Q — pYjI 
p 2 +0 -p) 2 


2 J2 


= <X 


because £[* 1 * 2 ] = 0 for independent variables, and 
fjjt j 2 j = fj JC2 j = cr ^ . The variance of the combined 
distribution then remains o 2 . Gaussian interpolation is 
therefore given by 


A w m,n = 


Pm.n Aw R{ 

^ nun ) 

i+i 

}~Pm.n ) 


A 

2 

Pm,rt 

+ i 

(l — Pm,n j 
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This combination is applicable for all three translational 
velocity components at an element location. 


The Dryden Model 

From reference 4 the low-altitude vertical turbulence 
scale length is given as a piecewise continuous function 
of altitude h 


The variance is defined by 


E[x]~ = x^f(x)dx = a 4 ' 


10 

/i < 10 ft 

h 

1 0 < * < 1 000 ft 

1000 

h>\000 ft 


for a single random variable. When we combine two 
independent random variables at some interior point, we 
must consider their proportion p. 

When a proportion p has the range (0 < p < 1 ), as it does 
in p m n above, a linear combination of two normally dis- 
tributed variables may be forced to have the same statisti- 
cal properties as the originals by normalizing the density 
functions as follows: 


/(*l) = 


f{ x 2 ) = 


7 2 J M 

Vp 2 +(1-p) 2 

/2 '~ p 2 /(*) 
Vp" + (!-p) 2 


and the in-plane components are given by the functions 


^14 


75.64 h< 10 ft 

hf~ v ] 2 1 0 < * < 1 000 ft 

1000 * > 1 000 ft 


where 

f uv =0.177 + 0.000823/? 

The horizontal turbulence RMS intensities are a function 
of the selected vertical turbulence intensity cr h . 

^ ,-0.4 

~ <* v — <J w f uv 

The results shown in this paper are for a w - 5 fVsec, 
h = 200 ft, and v H = 1 6.9 ft/sec ( 1 0 knots). 


As applied to variables with a zero mean value, this is a 
consequence of the arithmetic mean theorem given in 
reference 7, where 


Using Taylor’s frozen field hypothesis from reference 4, 
the Dryden form of the linear turbulence transfer func- 
tions are given by 
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For conventional aircraft (and our baseline configuration 
of the UH60A) the angular transfer functions are 
developed from partial differential equations relating 
these translational disturbances to the immersion rate of 
the vehicle into the turbulent field. These equations use a 
parameter b , which is the wingspan length of the aircraft. 


constant over each computer cycle At. The z-transforms 
of Laplace functions f£s) then become 


FM = z\ 




/«(*)! 


which may then be converted to difference equations for 
discrete implementation. This process is nicely explained 
in reference 8. 

Using this technique on the Dryden transfer functions 
produces the coefficients 

Yu = v uv^/^u 

fx=e -yu 

f2=0u('-f\h r Vn 

Yv = VuyAt/Lv 

gl = 2e~ Yv 


Ap(s) = 




s + 


4 b ) 


Aq(s) = 


4 b 


s + - 


TTVj 


4b 


83 = ][ I_e 7v +(S-\)e r 'Yv] 

84 =-\a v e ~ rv /4n][ ] ~ e ~ Yv +(V3~ I )y v j 

Yw - v uv ^/ ^vv 
h\ = 2e~Y' v 


“7*4v(s) 


,Y + — ^ 
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These angular disturbances are not used in the SORBET 
model because the geometry is considered in the points ot 
application of the translational turbulence excitations. The 
asymmetric turbulence velocities produce moments, 
which in turn produce angular activity. 

All of the baseline case transfer functions have been 
presented because their spectra are displayed on certain 
graphs herein for reference purposes. Also, this baseline 
model is actually available as an option in our UH60A 
simulation model, and was used by pilots tor comparison 
purposes. 

Discrete Implementation 

Continuous transfer functions involving random inputs 
are typically converted to discrete form using the zero- 
order hold formulation, where the input is assumed 


/? 2 = - e 

h3={°w/4y^ ][ 1 -<’ _r " +(V3-1 )e~ y »Yw] 

*4 = - [owe' 7 ” A [Y^ } [ I - e' 7 * + (' V3 - l))V ] 

These coefficients are used in the following six difference 
equations, each driven by an independent Gaussian noise 
source T] n with a zero mean value and a unity standard 
deviation. 

Au L (k) = f l Au L (k-\) + f 2 T 1] {k) 

Au R (k) = - 0 + f2 r l2( k ) 

Av L {k) = g\Av L (k - 1) + g 2 Av L {k - 2) 

+ 83*1 3( k ) + 84*1 3(k-l) 

Av R {k) = g\Av R (k - 1) + g 2 Av R {k - 2) 

+ g3n4(*) + S4'U(*- 1 ) 
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Aw L (k) = h\Aw L (k - 1) + h 2 Aw L (k - 2) 

+ / i 3 , ? 5 (*) + ^ 4 7 ? 5(* - 0 
Aw R (k) = h\Aw R (k - t) + h2Aw R (k ~2) 

+ h 3 ri 6 (k) + h 4 ii(,(k-\) 

These computationally efficient equations produce stable 
outputs for all stable inputs. 

For the body-fixed formulation these velocities are super- 
imposed at the center of gravity. In SORBET these 
velocities are created at the two onset points and then dis- 
tributed and superimposed at the various blade-element 
locations. They produce forces and moments because they 
modify the angle of attack of each element. 

Element Velocities 

Forces developed at individual blade-element locations 
are invariably summed to create total rotor forces and 
moments. In this section the spectral filtering conse- 
quences of the summation operation are illustrated by 
using the turbulence velocities themselves. Then, in the 
next few sections, this is further illustrated using the 
blade flapping spectra and the total vehicle responses. 
Although a blade variable may exhibit much harmonic 
content, the summation over the blade index effectively 
removes many of the harmonics (ref 9). 

Figure 2 displays typical time histories of vertical 
turbulence produced at the two onset points (Aw L ,Aw R ). 
These histories are the uncorrelated outputs of the 
Dryden vertical filters, excited by Gaussian noise with 
ct k = 5 ft/sec. 



The velocities of figure 2 are the origins of the turbulence 
velocities for each blade element. In response to these 
inputs, figure 3 is an example of the velocity of just one 
element (outboard, m = 5) of one blade. 


i5-r 



0 10 20 30 40 

Time - sec 

Figure 3. Element velocity. 


Because the velocities at an element are rotationally 
sampled, figure 3 displays considerable one-per-rev fre- 
quency content, bounded approximately by the input 
curves shown in figure 2. 

Each and every element experiences vertical velocities 
similar to that shown in figure 3. However, the summa- 
tion of all of the 20 elements for each time point selec- 
tively eliminates most of the high frequency content. 
Figure 4 displays this phenomenon. What remains after 
summation is essentially the average of the original onset 
velocities. 



Figure 2. Onset point velocities . 
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This phenomenon occurs independent of turbulence. In 
the next three sections this is further illustrated using 
both rotor and vehicle variables. 

Rotating Frame Blade Spectra 

The flight condition documented here is for a trimmed 
vehicle in 10-knot horizontal flight at an altitude of 200 ft. 
Under these conditions, when turbulence is not present, 
any given blade produces the flapping (p n ol ref. 1 ) time 
history as shown in figure 5. 



Time - sec 

Figure 5. Flapping without turbulence. 

The power spectral density (PSD) or autospectrum ol this 
signal is shown in figure 6. 



- 2-10 12 

10 10 10 10 10 

Frequency - Hz 


Figure 6 . Flapping PSD without turbulence. 

The significant frequencies at the blade level ot 
examination are the blade harmonics, beginning with the 
fundamental, at 4.297 Hz (27 rad/sec). These harmonics 
theoretically extend to infinity. Frequencies on this graph 
that are observed below the fundamental are, for the most 
part, aliased harmonics (originating beyond the Nyquist 
frequency). 


When the vertical turbulence {O w = 5 ft/sec) is introduced 
into the rotor system under the same conditions, the time 
history of a blade displays more frequency content, as 
indicated in figure 7. 



Time - sec 

Figure 7. Flapping with turbulence. 

The background spectrum of figure 6 is then elevated by 
a few orders of magnitude, as shown in figure 8. 
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Figure 8. Flapping PSD with turbulence. 

The significant periodic components in these spectra are 
created by the blades essentially tracking each other, 
independent of turbulence. Turbulence elevates the back- 
ground spectrum, but since the inputs are not uniform 
over the rotor disc it does not have a significant influence 
on the rotor harmonics. Notice that in figure 8, except tor 
the ftrst harmonic, the blade harmonics almost disappear 
into the turbulent background, and they are not amplified 
by the introduction of turbulence. This leads to the 
conclusion that at the blade level of consideration the 
introduction of turbulence does not produce significant 
cyclostationary phenomena. 
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The lagging spectrum is not shown here. At this flight 
condition it is about an order of magnitude smaller than 
the flapping spectrum. 

Fourier Coefficients 

Forces and moments generated by blade activity are 
transmitted to the airframe through summation functions. 
Given that each simulated blade has identical physical 
properties, such as mass, length, and relationship to adja- 
cent blades, the consequence of summation is to eliminate 
all integer harmonic multiples of the RPM that are not a 
multiple of the number of blades (ref. 9). This is illus- 
trated by considering the principal flapping Fourier 
coefficients. 

The principal flapping Fourier coefficients (main distor- 
tions of the rotor plane, in nonrotating frame) are given by 
the following summations of the individual flapping 
angles ( fi n of ref. 1 ). 

N 

Aof 

N 

A] F =_ ^ X ^ C0Sl// ' 2 

n = 1 
N 

B] f siny/n 

n = \ 

where Aq F is the coning (steady flapping) angle, A \ F is the 
longitudinal first harmonic, and B ]F is the lateral first 
harmonic. 

The PSDs of these Fourier coefficients are given in 
figures 9-11. The lower lines show the behavior without 
turbulence, and the upper lines show the behavior with 
vertical rotor turbulence. In comparing these figures with 
figures 6 and 8 it is seen that only blade harmonics that 
are a multiple of the number of blades survive the sum- 
mation operation, whereas the stochastic contributions are 
retained in the spectra of the rotor disc. 
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Figure 9. Aq F spectra. 
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Figure 10. A, F spectra. 
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Figure 11. B^ F spectra. 





Harmonics of N/rev seen in these figures are not caused 
by turbulence (see ref 9). They are a normal consequence 
of blade summation of periodic phenomena. Of course, 
interior harmonics will be observed if the blades are not 
tracking each other in their rotational environment. This 
occurs during transients, and is also a consequence of 
mismatched blades: 

The one-per-revolution (27 rad/sec) frequency 
content that can be seen in the flight-test data is 
caused by a mismatched rotor blade that was not 
tracking with the other three on the test aircraft; 
therefore, this frequency is not observed in the 
simulation results, (ref. 10) 


Baseline Traces 

The baseline spectra are the responses for the conven- 
tional implementation of turbulence, typically used for 
fixed-wing aircraft but also available in our UH60A 
model. In the baseline configuration the inputs are at the 
vehicle center of gravity. The transfer functions for this 
formulation were presented in the section named The 
Dryden Model. The baseline spectra also describe the 
outputs of SORBET’s right and left filters. 

The baseline configuration has been available for the 
UH60A model (and other helicopters) for a number of 
years, although pilots have invariably been critical of its 
performance during low-speed flight. This baseline 
configuration is nonetheless interesting for comparison 
purposes. In figure 12 the vertical velocity baseline 
spectra are shown, with parameter a = v/L w . 



For the flight condition of 10 knots at an altitude of 
200 ft, the parameter is about 0.04. 

At the 10-knot flight condition the baseline PSDs are 
presented in figures 13-18 as dotted lines. In the config- 
uration studied here only the vertical axis of SORBET 
had an input. However, all axes (solid lines) display 
outputs. 

N 




Figure 12. Baseline spectra. 
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Frequency - Hz 
Figure 15. Aw b , vertical input. 

It is useful to think of the body-fixed spectra (dotted 
lines) as representative of the velocity of the air mass , as 
experienced by a point traversing through it. If this air 
mass velocity is rigidly attached to the vehicle, as it is in 
the conventional formulation, then the vehicle itself 
experiences this velocity. If, however, it is attached to the 
velocity of a blade element, then it produces an alteration 
in that element’s angle of attack, which in turn con- 
tributes to the production of rotor-system forces and 
moments. 

The body-fixed baseline spectra could only be the same 
as the helicopter spectra if unity transfer functions 
existed from the rotor to the body. Indeed, the existence 
of nonunity transfer functions in this regard is obvious 
because rotorcraft responses differ from those of conven- 
tional aircraft. Nonetheless, these spectra are useful in 
determining just what the rotor-to-body effects are. Only 
“for high speed flight, where helicopters behave more like 
conventional aircraft” 2 could these dotted baseline curves 
possibly represent desirable responses. 

For completeness the angular PSDs are also shown in 
figures 16-18. All angular disturbances observed in 
SORBET were induced responses from the application of 
vertical turbulence. The baseline curves on these plots 
were obtained from the conventional transfer functions, 
using an arbitrary value of wingspan length b of 20 ft. 



Frequency - Hz 
Figure 16. Ap b , vertical input 


N 



Frequency - Hz 


Figure 17. Aq b , vertical input. 



2 Remark by Norman D. Ham of Massachusetts Institute of 
Technology during his class on Aerodynamics, Stability, and 
Control of Rotorcraft and Other VTOL Aircraft. Ames Research 
Center. 1978. 
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In figures 13-18 only the vertical component of turbu- 
lence into the rotor system was used, as was the case 
during pilot evaluations. Induced effects are seen to occur 
on all other axes. Pilots generally concurred that this 
system felt quite realistic, with the exception that some 
airspeed and lateral variations were lacking. This is illus- 
trated in figures 13 and 14. 

Rotor- Airframe Transfer Function 

In figure 15 the vertical vehicle response to vertical rotor 
turbulence was shown. Since the dotted line represents the 
spectral input throughout the rotor disc, clearly much 
energy is lost in the rotor-to-body transfer function. This 
transfer function was determined by emulating a single- 
input/single-output (SISO) system using the same 
statistical value for each element at each time point. 

The baseline vertical Dryden PSD as shown in 
figure 15 is plotted in figure 19 as a dashed line. Using 
identical random inputs at each element, the SISO PSD 
curve of figure 19 was created. Comparing the SISO 
curve to the input spectrum of the air mass (dashed line), 
there is considerable attenuation. Also shown in figure 19 
is the multiple-inputysingle-output (MISO) PSD curve, as 
replicated from figure 15. More attenuation occurs on this 
curve, when the inputs are not uniform (MISO) through- 
out the rotor disc. 



Frequency - Hz 

Figure 19. Aw b , same vertical inputs. 


The rotor-to-body transfer function (nondimensional) is 
also shown in figure 19, as computed from the SISO 
curve with respect to the baseline curve. The attenuation 
becomes more pronounced beyond about 4.3 Hz, which 
happens to be the rotor RPM. This phenomenon is a con- 
sequence of the blade summation process operating on 
signals that are identical for each element. 


When turbulence was permitted to vary over the rotor 
disc according to the algorithms, the slightly different 
results of figure 15 were produced, shown in figure 19 as 
the MISO curve. The transfer function for the MISO 
system has not been identified, although it must be sim- 
ilar to the SISO system given in figure 19 because of the 
vehicle response similarities shown by the SISO and 
MISO curves. 

Pilot Opinion 

From figure 19, the rotor system clearly attenuates high- 
frequency nonuniform inflow (as in blade-element 
turbulence) to a greater extent than it attenuates uniform 
inflow. Nonetheless, lower frequency distortion of the 
rotor disc caused by turbulence may be an important 
piloting cue. Distortions are an indication of activity in 
the stability augmentation system (SAS), and may influ- 
ence speed selection. Some of the simulation pilots stated 
that in real helicopter flight they sometimes “observed the 
rotor plane undergoing large responses to turbulence, 
while the airframe’s motion was relatively smooth.” In 
response to pilot requests for this visual cue, the tip-path 
plane will be included in the visual display for the next 
rotorcraft simulation on the VMS. 

To assure that simulation pilots had similar conditions in 
which to assess the model, they were given identical tasks 
of following recorded flightpaths. Our chief test pilot lor 
the simulation established these basic paths (without 
turbulence) by flying a set of runs over the terrain, and the 
spatial positions and orientations of the helicopter were 
recorded. Using a sophisticated playback feature of our 
heads-up display, other pilots were shown these paths by 
using both a graphical lead aircraft and a visually super- 
imposed stick-figure canyon (ref. 1 1). This is depicted in 
figure 20, where both the lead aircraft and the velocity 
vector (circle) are indicated. 



Figure 20. Playback canyon. 
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Pilots had the task of following the lead aircraft (at speed) 
and staying within the canyon. Without turbulence, this 
task was quite simple for all pilots. The task became more 
difficult when turbulence was added. 

All pilots found the SORBET model realistic, and their 
comments were used to improve the model during the 
simulation. 

Pilots generally agreed that the turbulence magnitude 
should be a function of velocity. Above 40 knots, the 
subjective values for light and moderate turbulence were 
selected as a w = 5 ft/sec and 8 ft/sec, respectively. Near 
hover, pilots selected standard deviations that were half 
these values. The reason for this is unknown, and should 
be investigated further. Although the frequency content is 
acceptable, the magnitude must be decreased in low- 
speed flight. For the body-fixed formulation in low-speed 
flight, pilots invariably criticized the baseline configura- 
tion. A typical comment was “It was like a washboard.” 

Pilots are clearly sensitive to the more nonlinear aspects 
of atmospheric turbulence. From reference 12, “turbulent 
flows are diffusive and intermittent” and, moreover, “the 
element of surprise that is characteristic of turbulence 
must be present.” Pilots were quick to point out that al- 
though the stochastic turbulence was realistic, it was too 
continuous and lacked occasional sharp gusts. Their 
comments reflected reference 13 quite closely: 

One of the most common grievances helicopter pilots 
have with simulated atmospheric turbulence is that it 
is too regular, and does not contain ‘patches’ of 
turbulence followed by periods of calm air. 

Due to the large number of pilot remarks in this area, we 
accommodated them on certain points. During the sim- 
ulation, a few additions to the SORBET model were 
made, as described below. The pilots agreed that these 
additions provided even more realism. 

Occasional vertical gusts were simulated by changing the 
magnitude of the mean vertical wind at random times. 
Gusts were programmed to arrive more frequently and 
more abruptly at high speeds (about every 3 sec) than at 
low speeds (every 12 sec). The resulting random vertical 
gusts were praised by the pilots as extremely character- 
istic of real turbulence. 

Moments of calm air were introduced by scaling the 
magnitude of the stochastic turbulence to zero for random 
time periods (between 2 and 6 sec). This was followed by 
smoothly applying random percentages of the full magni- 
tude, also for random time periods. This technique created 
the effect of patches of different levels of turbulence. 


With the addition of vertical gusts and patches of 
turbulence, the vertical turbulence was accepted by all the 
pilots as very realistic. One deficiency they found was in 
the reduced level of lateral, longitudinal, and yaw 
turbulence cues. Since the model did not directly influ- 
ence these axes, these comments were expected. As a 
partial solution, a stochastic input, perpendicular to the 
tail rotor, was added. While the pilots noticed the extra 
yaw variations, they generally agreed it was not a great 
improvement, and the lack of lateral inputs detracted 
from realism. 

One major pilot concern involved the way the turbulence 
was not correlated to the terrain. This was not within the 
scope of SORBET, but work is being done in this area 
(ref. 14). 

In summary, pilot comments were quite favorable for 
SORBET during the 2-week motion-based simulation 
study. Representative comments are listed below, where 
the italics are the authors’. 

1 . The chief NASA test pilot for this simulation, Munro 
Dearing, commented during low-speed flight that 
“this would be typical of light-to-moderate gusts in 
hover.” However, at 100 knots he noticed “there is a 
lack of heave.” These comments were made prior to 
the introduction of the embellishments described 
above. 

2. NASA test pilot Bill Hindson commented on 
SORBET once the occasional gusts were included 
that “the low frequency pitch and roll/*?// is also 
characteristic of turbulence,” and “it is a fairly 
realistic simulation.” 

3. NASA test pilot Tom Reynolds commented on 
SORBET when the occasional quiet periods were 
introduced that “this is much like what I expect is 
real,” and “overall, this feels pretty good, fairly 
realistic light-to-moderate turbulence .” 

4. NASA test pilot Gordon Hardy commented on the 
SORBET simulation that “ these are pretty typical 
excursions compared to real flight,” but also said “I’d 
expect more lateral motion.” 

5. Retired NASA test pilot Ron Gerdes commented 

“ there is nothing that I would change in the model” 

6. Navy pilot Tim Sestak commented “this is what 
turbulence feels like.” 

In-Plane Turbulence 

Although piloted runs were not made with the in-plane 
longitudinal and lateral turbulence options enabled, it has 
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been established that their inclusion would have added 
very little to the simulation. From reference 15, “Current 
articulated and semi-rigid rotors are insensitive to 
in-plane gusts.” Our own data support this statement. 

In fact, adding these in-plane components to SORBET 
produced changes that were too trivial to show here. For 
this reason we postulate that in-plane translational 
turbulence components should probably act at the vehicle 
center of gravity, just as in the conventional formulation. 
Of course, using these in-plane components also means 
that consideration should be given to correlated tail rotor 
responses. The tail rotor becomes much more significant 
when gusts are considered, rather than just stationary 
turbulence alone. 

In our next investigation we plan to include in-plane 
components for turbulence that will probably be applied 
at the vehicle center of gravity. This proposed system has 
also been investigated; it produces almost identical 
spectral responses as those shown in figures 13-18, 
except that the in-plane vehicle responses of figures 13 
and 14 then become superimposed with the dashed lines 
(as in the conventional aircraft case). This model should 
then include the in-plane disturbances that pilots expect. 

Conclusions 

Adding turbulence to the Black Hawk simulation added 
less than 10 percent to the required cycle time. Hence, 
blade-element implementations of turbulence are now 
feasible using modern computers. 

In-plane turbulence applied at the rotor elements produces 
minimal vehicle responses. This is therefore an unsatis- 
factory technique for modeling stochastic turbulence 
in the horizontal plane, and would probably produce 
erroneous results using gust models. Hence, in-plane dis- 
turbances should probably be handled more or less as in 
the conventional formulation, with the inclusion of the 
correlated influence on the tail rotor. 

The vertical component of turbulence as input at the blade 
radial stations produces significant differences in vehicle 
performance from that of conventional body-fixed tech- 
niques. The differences include considerable attenuation 
of higher frequencies in the vertical and all rotational 
axes. These phenomena elicit favorable pilot opinion, and 
probably identify the major difference between helicopter 
and conventional aircraft responses. 

In establishing a uniform level of turbulence over a flight 
regime our test pilots preferred a standard deviation that 
varied linearly with velocity, roughly doubling from 
hover to 40 knots, where it became constant. This is an 
unknown phenomenon that deserves further investigation. 


Also, criteria for implementing patches of different levels 
of turbulence should be established, because they 
correspond to pilot experience. 
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